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- Mass diffusion, including self-diffusion, mutual diffusion, tracer diffusion, etc. 
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Abstract - The paper presents the velocity autocorrelation spectrum of bulk water measured by 
a new technique of NMR modulated gradient spin echo method. This technique is unprecedented 
for the spectrum measurement in the frequency interval between a few Hz to about 100 kHz with 
respect to directness and clarity of results and shows that a simple model of Brownian self-diffusion 
is not applicable to describe the diffusion dynamics of water molecules. The observed temperature 
dependant spectra of water show the existence of a slow chain-like dynamics in water, which we 
explain by coupling of diffusing molecule to broken bonds in the hydrogen bond network. 



Water properties are essential to many physical, chemical, biological, and geological 
processes. Although a great deal is known about thermodynamic properties of bulk water, 
the complex forces that govern its molecular arrangements and dynamics obscure a detailed 
picture of instantaneous local structure of water and its evolution. Understanding these 
phenomena has been one of the most important scientific challenges of the past 100 years. 

In water, intcrmolccular hydrogen bonds lead to a highly associative structure that has 
attracted the interest of scientists because it exists in the temperature range associated 
with the existence of life. According to current knowledge, the structure of liquid water 
is a disordered network of molecules connected by hydrogen bonds, which behaves on time 
scales of 10~ 12 to 10~ 9 sec as a gel. The structure fluctuates and reorganizes, due to 
rotations and other thermal motion, causing the unceasing breaking of individual hydrogen 
bonds and their reformation into new configurations, whose extent and influence depend on 
temperature and pressure. It is thought that the hydrogen-bond network does not remain 
intact for long enough to be detected directly as an observable entity in ordinary bulk water 
at normal pressures. 

Insights into structural properties of water have come from x-ray and neutron-scattering 
experiments, which lack dynamical information; insights into its dynamics have been gained 
from ultrafast time- resolved experiments, which have lacked structural details [TJH]. The 
most detailed understanding of water properties derives from molecular dynamics simula- 
tions, which provide an atom-by-atom perspective on how hydrogen bonding changes with 
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time. Although substantial efforts have been made to develop models for molecular dy- 
namics simulations [3] that show some degree of success in reproducing the experimentally 
observed radial distribution function pQ, the models have difficulty addressing the unusual 
nature of long-range ordering that strongly affects the properties of water and dynamics of 
ordering has never been properly tested against the experiment. 

The traditional experimental techniques are yielding diminishing returns, owing to the 
collective efforts of many generations of researchers. Instead, new results should be expected 
with techniques unrelated to these mainstream research efforts. For example, the recent 
dielectric measurements show an anomalously slow relaxation process in water [5], which 
is believed to be a kind of collective motion of hydrogen-bonded structures. It gave an 
impetus to apply other nontraditional techniques for the study of long range dynamics in 
water, among which the measurement of translation dynamics by a new NMR technique 
can provide valuable information. 

Many of the equilibrium and dynamic properties of liquid water show anomalous behavior 
that is changing with temperature and pressure. Example is translational self-diffusion, 
which is a key property for a long range molecular dynamics. Its measurements by tracer 
methods [6], NMR techniques [TIE] and indirectly by non-elastic neutron scattering |9] 
reveal anomalously strong temperature dependence of self-diffusion coefficient that can be 
attributed to the intermolecular hydrogen bonding in liquid water [TU1[TT]. However, the 
diffusion coefficients calculated for many theoretical models of liquid water give mostly 
larger values than experimentally observed |12j . 

NMR pulscd-gradient-spin-echo NMR method (PGSE) plays a prominent role in the 
measurements of the molecular self-diffusion. PGSE determines the molecular displacement 
in a time interval between two pulses of applied magnetic field gradient via the spin echo 
attenuation [13j . However, the results of measurements in water at different temperatures 
and pressures differ among authors [51t51ll4j. on a scale exceeding the experimental uncer- 
tainty. Presumably the uncertainties in the gradient calibration and convection flows due to 
temperature variation in the sample are the cause. The measurements with the new NMR 
technique presented here show that the scattering might be a result of the time-dependent 
self-diffusion coefficient, i.e. non-Brownian diffusion of water molecules [15j that has been 
overlooked in the PGSE experiments. Nevertheless, PGSE measurements brought forth the 
poor validity of the hydrodynamic model of water [5]. Namely, the obtained self-diffusion 
coefficient D expressed in terms of the Stokcs-Einstcin-Dcbyc formula of particle diffusion, 
D = , shows a dramatic change in the effective radius R of the diffusing species in the 
temperature interval from 300 to about 600 K. This necessities a new molecular description 
of self-diffusion in water that takes intermolecular hydrogen bonding into account. 

A step in this direction comes from a new viewpoint on the labelling of moving spin- 
bearing molecules by a sequence of magnetic field gradient and radio-frequency (RF) pulses 
in NMR. It turns out as a natural relation of the accumulated time dependent spin phase 
to the velocity autocorrelation spectrum (VAS) of molecular translation motion |16| . which 
is defined as [T7] 



where Av z (t) = v z (t) — (v z (t)) is the component of the velocity fluctuation in the direction 
of the applied gradient (z-axis in our case). Repetitive sequence of RF and magnetic field 
gradient pulses results in a periodic spin phase modulation. The frequency spectrum of 
this modulation has a narrow peak, which can be used to sample VAS by changing the 
modulation frequency [18] . This technique of modulated gradient spin echo (MGSE) acts as 
a filter with an adjustable detection window into the frequency domain of spin translation 
dynamics because it samples only molecular motion that is in pace with the spin phase 
modulation. The method is unique in a clear discrimination of fast diffusion-like motion 
from slow collective one, e.g. flow or convection in fluids. Therefore, it does not require a 
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special temperature stabilization to prevent convection flows within the sample, which is a 
big obstacle in the PGSE measurement. VAS of ordinary liquids decreases at frequencies 
higher than the reciprocal correlation time of velocity fluctuations being in the pico-second 
range. Its constant value at low frequencies determines the Brownian character of molecular 
self-diffusion [THIHO] . MGSE measurements of liquid reveal whether character of molecular 
self-diffusion is Brownian or not and cast a new light on the structural dynamics. 

The pulsed-gradient version of the MGSE method has already been used to study flow 
and self-diffusion in porous media [2"lTl2~i] . but the inductance of gradient coils limits the 
switching time of gradient pulses to the modulation frequencies below 1 kHz. For this 
reason, the technique does not attract a lot of attention. For instance, in porous media, it 
cannot explore the characteristic of VAS typical for restricted diffusion in pores smaller than 
1 /xm. Recently, we demonstrated [25] that a standard Carr-Purcell-Meiboom-Gill sequence 
of 7T-RF pulses [55] applied simultaneously with a constant magnetic field gradient acts as 
a MGSE method. This technique overcomes the high-frequency deficiency of the pulsed 
MGSE technique, because the absence of gradient-switching induced eddies enables VAS 
measurement up to the frequencies of about 100 kHz. In the NMR receiving coil aligned 
along x-axis, spins induce signal 

E(r) - j t Y,Trp,(r)I xi , (2) 

i 

where the dynamics of a single spin operator I X i in the magnetic field, and the molecular 
motion govern the evolution of the density matrix 

Pi(r) = Ui(t) Po U+(t) = U im (t)U is (t) Po U+(W? m (t). (3) 

By additional factoring of the spin evolution operator f7j s (i), we detach the spin interac- 
tion with the main magnetic field Ui (t), the interaction with RF field UmF(t), and the 
interaction with gradient field G, [THIHI] governed by the operator 

Ui G {r) = Texp[ijG- Ti(t){I zi cos w a(t) + 

Jo 

+ I X i sin it a (t)}dt}. (4) 

Here T is the time ordering operator, r^(i) is the fluctuating location of the i-th spin, and 
the factor a(t) alternates between and ±1 with the rate of spin phase inversion by 7T-RF 
pulses [16| . The second term of the exponent includes the resonance offset effects, which 
have to be reduced as much as possible. Proper phase cycling of 7T-RF pulses suppresses the 
second term to an extent that the combination of RF pulses and constant magnetic field 
gradient behaves as an effective gradient with the sign alternating with the rate of applied 
7T-RF pulses, when the amplitude of RF field is larger than the difference of the non-uniform 
magnetic field across the sample. This condition implies that the gradient selected slice 
must not be smaller than the sample. The gradient spectrum, 



G(w, t) = G f cos(7ra(t))e ia,i dt , 
Jo 



has a peak at u rn = n/T m , where T m is the time interval between the 7T-RF pulses. The 
resulting spin echo attenuation [2"5] : 



8 7 2 G 2 r^(^ m ) 



(5) 
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Fig. 1: Measured velocity correlation spectra of water at 21° C (blue) and at 4° C (red) are 
fitted with the spectra calculated from the Langevin equations of harmonically coupled Brownian 
particles. The spectrum of ethanol at 21° C (green) is given to emphasize the difference between 
water and liquid with a weaker hydrogen bonding. 

is proportional to VAS at the modulation frequency 1/T m . By changing T m we can sam- 
ple D z (uj), from which we can obtain the mean-squared displacement of the spin bearing 
particles using the relation ([z(t) — z(0)] 2 ) = ^ J °° D z (uj)(1 — cos(ujt))/uj 2 duj. 

Our experiments were carried out on a TecMag NMR spectrometer with a 2.35 T hori- 
zontal bore superconductive magnet equipped with micro-imaging accessories and Maxwell 
gradient coils. Samples were sealed in a small glass tube of 0.7 cm diameter and 2 cm 
length. With this equipment, we have demonstrated the use of the new MGSE technique 
by measuring VAS of restricted diffusion in porous medium with pores smaller than 0.1 ^m 
[2"5] and by the first observation of VAS of fluidized granular system . 

In a study of self-diffusion in gel systems by this technique, water was meant to be used 
as a calibrator at first. However, much to our surprise VAS of water exhibited a very distinct 
low-frequency decrease in the range below 2 kHz as shown in fig. [T] Systematic error was 
ruled out by identical measurements of other liquids: nitrobenzene, toluene and ethanol, 
which all yielded a flat VAS in this low frequency range (fig. [J). Possibility of water heating 
by a high duty cycle of RF pulse train was ruled out by performing experiments with the 
same modulation frequencies but 40 ms and 400 ms long CPMG RF trains. Both measure- 
ments give identical frequency dependencies. The expected result of MGSE measurements 
is also a strong temperature dependence of water VAS in the temperature interval of our 
measurements, from 4°C to 40°C, which at low frequencies converge within experimental 
error towards the results of Mills. Mills measurements by the diaphragm technique are 
commonly considered as the most accurate measurements of long time water self-diffusion. 

The observed VAS of water has features similar to those of restricted self-diffusion in 
porous media [HJ[25] or of polymer self-diffusion [3D], which indicates interaction of the 
diffusing molecule with surrounding and thus constraining its displacements. As shown 
in the following, a simple model of a Brownian particle interacting with local attractive 
centers [3T] fits well the observed low frequency decay of VAS. Therefore, we assume that an 
association of the diffusing molecule with incessantly breaking bonds of the hydrogen bond 
network in water is responsible for its constrained motion. This also explains a stronger low 
frequency decrease of VAS at lower temperature, fig. [TJ as the structure of hydrogen bonds 
becomes more stable. 

This reasoning follows the lines of a widely accepted view of tetrahedral arrangements 
of water molecules ordered by hydrogen bonds, where the mechanism of fast switching be- 
tween the hydrogen-bonded partners remains unclear. Recent combined studies by 2D IR 
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spectroscopy and molecular dynamics simulations [32j[33] revealed that a typical configura- 
tion has many instances of molecules with apparently broken hydrogen bonds. Simulation 
studies |34j showed that although the majority of molecules in liquid water participate in 
four hydrogen bonds, many of them participate also in two, three or five hydrogen bonds at 
a given time. Recent observation of slow dynamics in the frequency range below 1 kHz by 
the dielectric measurement of water [5] indicated a possibility that the distorted hydrogen- 
bond structure behaves like a chain in polymers. Although translational molecular motion 
itself does not lead to a significant bond breaking or bond restoring, the results of MGSE 
measurement give a new insight into interactions of molecule diffusing within the hydrogen 
bond network in water. 

To interpret the distorted VAS in water, we use a simple model of Brownian particle 
harmonically coupled to attractive centers [311 [35] i.e., to broken bonds of the hydrogen 
bond network. We assume that instability of the network introduces the fluctuation of 
the attractive centers, which is independent of the forces acting on the diffusing molecule. 
As the relaxation times of the forces associated with the motion of libration are short in 
comparison to the times of the translational motion, the Markovian approximation of the 
stochastic process permits to describe the molecular displacements along the applied gradient 
with coupled Langcvin equations without the oscillating terms: 

-^ + K{Z-z) = g{t). (6) 

Here k is the coupling constant between the diffusing molecule, at the coordinate z, to 
the attractive point in the hydrogen bond network located at Z. £i and £2 are friction 
coefficients while f(t) and g(t) are stochastic forces acting on the diffusing molecule and on 
the attractive point respectively. Calculation gives the low frequency part of VAS as 

1 + ^& 2 T C 2 

DH = -D ° 2 2 • ( ? ) 
1 + uj z t£ 

Only three parameters characterize the spectrum: the characteristic time of molecular cou- 
pling to attractive centers r c = /^^n , the high-frequency plateau is the diffusion rate 

resulting from short almost unconstrained molecular jumps D(oo) = </(t l 2>Tc , and the dif- 

22 

fusion rate at long displacements D(0) = ^^^^^i ^ Tc , is affected by the network inter- 
action. For fixed centers, the value of D(0) is zero but, in our model, the random breaking 
and restoring of bonds results in a finite value of D(0) and is a measure of network stability. 

By assuming the probability that the diffusing molecule can form multiple bonds with 
the free nodes on the network during its travel through the medium, we have obtained the 
model spectrum that fits well to the experimental data as shown in fig. [1] At 4°C, the best 
fit is obtained if the molecule is coupled 40% of the time with a force four times weaker then 
the coupling force for the remaining 60% of the time. At 20°C the best fit is obtained if 
the coupling is mostly of the weaker type. Fit with intermediate values of coupling deviates 
from the experimental data. We attribute the four-to-one ratio of the forces to the molecular 
association with one hydrogen bond for the weaker type and with four hydrogen bonds for 
the stronger one. Thus, effective interaction of diffusing molecule at room temperature is 
one-hydrogen bond, while at 4°C it is the combination of one and four hydrogen-bonds 
in almost equal proportion. VAS at 40° C (not shown here) is almost identical to that 
at 20°C, but with the plateau lifted to 2.9 x 10~ 9 to 2 /s, what could mean that one-bond 
interaction is dominant at higher temperatures, while the tetrahedral arrangements prevail 
near the freezing point. Plots of measured spectra as a function of ^Juj show a linear 
dependence at very low frequencies, which might indicate a chain-like motion similar to 
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that in polymers [36], when only one hydrogen bond is involved. According to this model, 
a larger difference D(oo) — D(0) at 4 ° C compared to those at higher temperatures could 
indicate a higher network stability at low temperatures. 

Evidently, the results of MGSE measurement of the bulk water open some interesting 
questions about the impact of the network-like structure on the molecular self-diffusion. 
A simple theoretical explanation for VAS of water is provided, however, a more detailed 
theoretical analysis still needs to be done. 

* * * 
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